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We propose and discuss a scheme for robust and efficient 
generation of many-particle entanglement in an ensemble of 
Rydberg atoms with resonant dipole-dipole interactions. It is 
shown that in the limit of complete dipole blocking, the sys- 
tem is isomorphic to a multimode Jaynes-Cummings model. 
While dark-state population transfer is not capable of creat- 
ing entanglement, other adiabatic processes are identified that 
lead to complex, maximally entangled states, such as the N- 
particle analog of the GHZ state in a few steps. The process 
is robust, works for even and odd particle numbers and the 
characteristic time for entanglement generation scales with 
N a , with a being less than unity. 



Entanglement is one of the most distinct quantum 
features of many-particle systems and has only recently 
started to be studied in a more systematic way. It pro- 
vides strong tests of quantum nonlocality pi and is at 
the heart of quantum information science with numerous 
applications [0-||. One of the open practical problems 
is to identify mechanisms for its robust and controlled 
generation. Recently Maimer and S0rensen suggested 
a scheme to create the TV-particle analog of the GHZ 
state in an ion-trap without the need for a precise con- 
trol over the collective vibrational modes of the ions || . 
Due to Kramers degeneracy M in the underlying nonlin- 
ear Hamiltonian different unitary operations needed to be 
applied for even and odd number of particles. The opti- 
mum interaction time T scales linearly with the number 
of particles and extreme fine tuning of T is required. As a 
consequence this method is highly sensitive to variations 
of external and internal parameters and in the presence 
of decoherence the success probability decreases expo- 
nentially with N . 

In the present paper we discuss a robust and efficient 
method to create complex entangled states like the N- 
particle analog of the GHZ state m 
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in a 3-step adiabatic process in a total interaction time 
which scales substantially less than linear in the number 
of particles. The underlying interaction is the resonant 



dipole-dipole interaction and the associated blockade ef- 
fect in an ensemble of frozen Rydberg atoms studied in 
It was shown in |j] that the dipole-blockade can 
be used to generate any (symmetric) entangled many- 
particle state by applying a specific sequence of resonant 
pulses Jic]] . A substantial drawback is however the need 
of a large number of pulses (which scales linear in N) 
with well defined pulse area. Hence the method is also 
highly sensitive to variations of external and internal pa- 
rameters and due to the intermediate excitation of de- 
caying Rydberg levels the success probability decreases 
exponentially with N. 

Following the proposal of Ref. || let us consider an en- 
semble of N atoms with two lower levels \a) and \b) both 
coupled to a Rydberg state \r) by coherent laser fields 
with (real) Rabi- frequencies f2i(t) and ri 2 (t) respectively. 
Let us further assume that there are two additional Ry- 
dberg levels above and below |r) with equal energy split- 
ting. In such a configuration there is a resonant energy 
transfer between two atoms in Rydberg levels, leading to 
a symmetric splitting of the doubly-excited states. If the 
minimum splitting, given by the atoms of largest separa- 
tion, exceeds the natural linewidth, resonant laser exci- 
tation into doubly- or higher excited states is suppressed 
(dipole-blockade). In the absence of this dipole-blocking, 
the Hamiltonian is linear in the total spin of the atoms 
(SU(2) symmetry) and it is not possible to create entan- 
glement. 

On the other hand for perfect dipole-blockade there is 
never more than one excitation in the Rydberg levels. 
In this limit the effect of the dipole-dipole interaction 
can easily be modeled by treating atoms in the Rydberg 
state as fermions (a, <J + ), while representing atoms in 
levels | a) and \b) by bosons (a, a + ), (6, b + ). The presence 
of the fermionic component breaks the SU(2) symmetry 
and the interaction is no longer linear in the total spin 
but can be described by a multi-mode Jaynes Cummings 
Hamiltonian 

H(t) = A ia +a + A 2 b+b + (pid + fl 2 bja + + h.c. (2) 

The detunings Aj have to be much smaller than the min- 
imum splitting of the doubly-excited manifold for the 
blockade- limit to hold. 

The isomorphism to the multi-mode Jaynes Cummings 
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model has a number of interesting consequences. First of 
all it simplifies the analysis by allowing to employ angular 
momentum techniques. Secondly many known features 
of the Jaynes-Cummings dynamics, such as decay and re- 
vivals of oscillations squeezed-state generation, and 
quantum state transfer between different modes |l3| ] can 
be anticipated in the dipole-blocking system. 

The blockade of double and higher excitations results 
in a chainwize coupling between symmetric collective 
states as shown in Fig. 1. This coupling with an odd 
total number of levels suggests the application of dark- 
state Raman adiabatic passage techniques [Q. To ana- 
lyze adiabatic passage in such a system it is convenient 
to introduce dark- and bright-state boson operators 

D = a cos 6 — bsin9, B — a sin 6* + 6 cos 6* , (3) 

with tan# = f2i/£7 2 . In terms of these variables the 
interaction Hamiltonian reads 
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+ n (B<j+ + B+ a), 

with Slo = \An(*) + W- The first two terms are the 
free energy of the atoms in the dark and bright states 
and the third term describes the coupling between dark 
and bright states. The last term shows that only the 
bright-state component is coupled to the Rydberg levels. 



population, the presence of dipole-dipole interactions and 
the resulting dipole blockade are irrelevant. Neverthe- 
less, as will be shown in the following, adiabatic tech- 
niques can be used to create entanglement if other than 
the zero-eigenvalue state are involved. 

To this end we consider here a situation opposite to the 
two-photon resonance when A = Ai = — A 2 . We first 
discuss the case of a substantial delay between the two 
pulses Hi and fi 2 such that the coupling between the dark 
and bright components, which is proportional to A sin 29 
is negligible. In this approximation the Hamiltonian (^) 
can be expressed in the simple form 



H = — a z cos 29 
2 
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where the irrelevant constant term A r A/2cos20 was 
omitted. The corresponding Schrodinger equation can 
be solved analytically in the adiabatic limit, i.e. when 
the mixing angle 9{t) changes sufficiently slowly in time. 

To obtain a convenient closed form of the solution we 
introduce angular momentum operators J\ = (a + B + 
a-B^)/2y/M,J 2 = i(a+B - a~B^)/2VM, and J 3 = 
a z /2, where M = B^B + a + a is the constant parti- 
cle number in the bright-state-Rydberg manifold. In 
terms of these operators the Hamiltonian reads H = 
Oe- i,3j2 J 3 e^ j2 , where U = V^oW + A 2 cos 2 20, and 
tan/3(i) = Q (t)/Acos29. The corresponding unitary 
evolution operator then reads 



U (t) = c - imj2 exp 



i,J 3 / n(t')dt' 



J/3(-co)J 2 



(0) 
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FIG. 1. Coupling scheme of collective iV-atom states in 
limit of dipole blockade, shown here for Ai = A2 = 0. In- 
dividual atoms have two lower states \a) and \b) coupled to 
Rydberg state |r) with Rabi-frequencies Qi and Q.2 respec- 



tively. 



z h m ) denotes symmetric superposition of N — m 



atoms in state \a) and m atoms in state |6). 



Under conditions of two-photon resonance, i.e. Ai = 
A 2 the dark-state subspace decouples from the remaining 
system. Its dynamics has however again SU(2) symmetry 
and factorized states remain factorized. Hence dark-state 
adiabatic transfer is not suitable for entanglement gen- 
eration. This result can easily be understood physically. 
Since the dark state does not contain any excited-state 



Let us now consider the case of all N atoms being 
initially in \a). If an intuitive pulse sequence is applied, 
i.e. if ill is switched on and off before f2 2 one has (3 : 
— > 7r and the systems starts from a bright state 
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where |a w ) denotes the collective state with all atoms 
being in level \a). The unitary evolution operator then 
reads 



W = -2iJ 2 exp 



+00 

-U 3 fdtp(t) 



(8) 



One-time application of W generates a symmetric collec- 
tive state containing a single Rydberg excitation and all 
other atoms are in \b): 
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corresponding to 



Applying W twice generates the W-state of Ref. fL5| 



| a r) and hence its amplitude undergoes Raman adi- 
abatic passage to |b Ar_1 r) through a chain of successive 
(10) T^-type transitions. Since the fields are not in iV-photon 
resonance, it is essential that If2 m | 3> |A|. 



(11) 



On the other hand starting from an initial state with all 
atoms in \b) corresponds to a pulse sequence in counter- 
intuitive order and leads to the transfer 
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Iterative applications of the same operator W allows 
one to reach any state in the 27V+1 dimensional manifold 
of symmetric many-particle excitations with at most one 
Rydberg atom. The VT-operation is based on adiabatic 
evolution and is robust against variations of parameters 
as long as the condition 
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is fulfilled, with 7 being the decay rate of the Rydberg 
levels. 

Although the application of W leads to an entangled 
state whose creation would require many 7r-pulses, O(N) 
steps are needed for the generation of complex states like 
the iV-particle analog of the GHZ state (Q) . We will now 
show that a small modification of the W operation can 
achieve this goal in very few steps and independent on 
the number of particles. 

For this we assume that the system is initially in an 
equal superposition of atoms being in \a) and the sym- 
metric state containing a single Rydberg excitation. 



a N ) + la^-V 
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|$o) can easily be created out of \^o) in a robust way 
e.g. by sweeping the frequency of f^i through resonance 
(rapid adiabatic passage) ||16). We now apply the W 
operation discussed above, however with a smaller time 
delay between the two pulses. In this case the dark- 
bright state coupling in the Hamiltonian (|]) proportional 
to sin 29 needs to be taken into account. Furthermore it 
is assumed that |fio| ^ Under these conditions the 
Schrodinger equation can no longer be solved analyti- 
cally. However numerically evaluating the equations (for 
N up to 20), we found the behavior shown in Fig. 2. 

The mechanism can qualitatively be understood as fol- 
lows: Due to the non-vanishing detuning A and the cho- 
sen intuitive pulse order, the state amplitude in \a N ) un- 
dergoes an adiabatic return process Jl6[ and ends up in 
the same state as it started from. At the same time 
the chosen pulse order is counter-intuitive for the state 
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FIG. 2. Temporal evolution of population in la^) and 
I b> 1 r) form initial state f° r N = 5. The laser pulses 
are Gaussian fli t 2 (t) = Cl m exp[— (t ± r) 2 /T 2 ], the delay is 
t = 0.5T the pulse area f2 m T = 125, and AT = 50. 



The amplitude of the state vector in la^) undergoes 
some rapid oscillations but returns to the same state for 



t 



00. At the same time the amplitude in |a 



JV-l, 



is transferred to |b r). Applying in a third step the 



inverse of W, eqs.(|ifj|) and (12), eventually leads to the 
iV-particle GHZ state (|l|). This corresponds to the overall 
3-step adiabatic process: 
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JV-l. 



The transfer is in all parts robust. It does not depend 
on the exact pulse form of Cl\ and O2, nor does it re- 
quire an extreme control of the delay time r or the pulse 
length T. Furthermore the mechanism works for even 
and odd numbers of atoms in the same way. In Fig. 3 
we have shown the dependence of the final population in 
the states la^) and |b iV ' _1 r) for N = 5 as function of 
pulse delay t and pulse area £l m T . It can be seen that 
the mechanism is robust against small variation of the 
delay time and - above some critical limit - the pulse 
area. It should be mentioned that for very large values 
of the pulse area the populations decrease again, since 
then the term A sin 29 in is negligible and there is a 



transfer -^(\ 



\a N - l r 



l&^-V)). 

Since all processes are adiabatic and Cl m is limited by 
the dipole-blockade condition, the question arises how the 
time T for generating the GHZ state scales with the num- 
ber of particles. From our numerical calculations, shown 
in Fig. 4, we find T ~ N a with a < 1 and decreasing with 
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N. The numerical calculations for N = 3 ... 16 indicate 
a — > 2/3. Thus in the presence of decay, the success 
probability decreases less than exponential with N. 
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FIG. 3. Final population of states a^) and |b Ar_1 r) for 
conditions of Fig. 2. top: as function of delay time r for 
Q m T = 120, bottom: as function of f2 m T for r = 0.5T. 
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FIG. 4. Minimum pulse area f2 m T m i n to create GHZ state 
as function of particle number N. Dots represent values from 
numerical solution of iV-particle Schrodinger equation for op- 
timized r. 



In conclusion, we have proposed an efficient and ro- 
bust method to generate complex entanglement struc- 
tures, such as the TV-particle GHZ state in a many- 
particle system with resonant dipole-dipole interactions. 
The method is robust against variations of parameters 
since for all steps adiabatic transfer processes are used. 
Although dark-state adiabatic passage is not suitable for 



entanglement generation, as it does not involve popula- 
tion of the interacting Rydberg levels, other adiabatic 
processes are identified that allow e.g. for the genera- 
tion of the iV-particle GHZ state (|l|) in three steps. The 
suggested method avoids the problem associated with 
Kramers degeneracy and thus works for even and odd 
number of particles. Exact knowledge of the number 
of particles is not required, making the method robust 
against number fluctuations. As opposed to the proposal 
of ref. H no extreme fine tuning of the interaction time 
is needed and the minimum interaction time scales only 
less than linear with the number of particles. Finally it 
should be mentioned that similar ideas can be applied to 
other many-particle systems, e.g. ions in a trap. 
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